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Abstract. 

The Canadian Climate Centre (CCC) GCM has been 
modified to allow its use for studies in atmospheric chem- 
istry. The initial experiments reported here have been run 
to test and allow sensitivity studies of the new transport 
module. The impact of different types of parameterization 
for the convective mixing have been studied based on the 
large scale evolution of Rn 222 and Pb 210 . Preliminary re- 
sults have shown that the use of a scheme, which mixes 
unstable columns over a very short time scale, produces a 
global distribution of lead that agrees in some aspects with 
observations. The local impact of different mixing schemes 
on a short lived tracer like the radon is very important. 


1. Introduction. 

Recently the modeling of chemistry has moved to the global 
scale, permitting the assessment of the multitude of chem- 
ical reactions within a global atmospheric model. The 
model described herein, the AES and York University Chem- 
istry and General Circulation model, is a version of the 
Canadian Climate Centre General Circulation Model, (GCM). 
The model used is the operational version of the CCC 
spectral GCM (McFarlaneet al, 1992). Initial experiments 
were aimed at examining the use of the spectral method for 
tracer transport and assessing the transport module. This 
module includes, multi-tracer advection, convective and 
boundary layer mixing, wet deposition, vertical and hori- 
zontal diffusion. This work describes some experiments on 
radon and lead tracers chosen to assess the part played by 
convective mixing in the transport of species within the 
GCM. 

The current results are for T32 resolution and the model 
time step is 20 minutes. At this resolution the model re- 
solves synoptic and planetary scale motions, however sub- 
grid meteorology which impacts on these scales are absent. 

We have chosen Rn 222 as being representative of species 
with short chemical lifetimes such as isoprene and propane, 
and Pb 210 as being representative of a longer lived species, 


e.g. CO and 0 3 . Very simplistically we are equating Rn 222 
with a precursor chemical like NO x or hydrocarbons which 
have short lifetimes, while the lead is equated to the longer 
lived chemical products, whose presence in the upper tro- 
posphere can only be achieved if its precursors are trans- 
ported rapidly from their surface source regions. In ad- 
dition, radon tests the impact of a simple surface source 
while the Pb 210 production is distributed throughout the 
troposphere. 

2. The Radon Experiment. 

As a first test of the model and transport module, a num- 
ber of paired Rn 222 and Pb 210 like tracers were used. This 
allowed several convective mixing schemes to be considered 
and a general assessment to be made of the overall use of 
the spectral GCM as an online tracer advection model. A 
test of the mixing, rainout and spectral advection plus a 
first check on mass conservation and other numerical prob- 
lems associated with sharp tracer gradients modeled in a 
spectral model was thus possible. 

Initially, radon was distributed uniformly in the hor- 
izontal, with a fixed vertical profile, decreasing exponen- 
tially with height. The ground level value of 100 pCi/m 3 
STP was chosen based upon typical observed surface radon 
measurements. Lead was set as a low background which 
would be dominated rapidly by the radon decay source. 
Due to its short lifetime, the radon rapidly reached a global 
balance between the model source and decay within the 
first month and was largely unaffected by the initial state 
thereafter. For the lead, a global mass equilibrium state 
was not expected wit hin the first several months due to its 
longer tropospheric lifetime. Since the radon surface flux 
is confined to the continent, the sharp boundaries and the 
hemispheric gradient, due to the source distribution, pro- 
vide test problems for the model transport and analysis. 
The radon source was modeled by assuming a constant 
surface tlux of 1 atom/sec/cm 2 from all non-ice covered 
land on the globe. The radon sink due to radioactive de- 
cay was obtained using a tendency based on the tracers 
5.5 day e- folding time. Thus excluding the transport, the 
budget, lor Rn 222 is described 
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in equation 1 , 


^I = _* L + F 
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where \ r = Rn 222 tracer mass mixing ratio, r= radon e- 
folding time, and F g = surface flux of radon. The Pb 210 was 
assumed as the only decay product of the radon so that 
the sink of radon in molecules was the same as the source 
of Pb 210 in molecules. Its sink was modeled using a simple 
rainout scheme, proportional to the local column rainfall 
amount and the local tracer mixing ratio (Mahlman and 
Moxim, 1978). Thus the budget for Pb 210 is; 
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where \ p = Pb 210 tracer mass mixing ratio, A/ p = molecular 
mass of Pb 210 , A/ r = molecular mass of Rn 222 and t to = 
rainout e-folding lifetime as defined below as; 

T U! (k'). < P > 

T.„ = ? (.5) 

where r w = wet deposition lifetime (as a function of height, 
k,), P = local average precipitation amount (mms -1 ), 
and < P > — the global average precipitation amount 
(nun ; 1 ). 

Convective mixing of tracer from the boundary layer 
into the free troposphere has been identified as an impor- 
tant component in a complete transport description of the 
troposphere. Tracer mixing schemes that have been used 
for radon include those of .Jacob and Prather (1990), l e- 
ichter and (Tutzen (1990) and Brost and (’lialfield (1989). 
These schemes vary in complexity with a balance being 
struck between detailed model description and expense. 
However duo to the large scale of the (KM s horizontal 
grid, parameterizations are, at best, highly approximate. 
Such a position poses a serious difficulty in global model- 
ing and requires careful assessment and observational basis 
before interpretation can be considered valid. 

It was decided to test, t wo schemes in parallel using fun- 
damentally different transport, assumptions. The two basic 
schemes used art;; a) a direct relocation scheme and h) a 
diffusive scheme. The direct relocation scheme is based 
on the scheme described in Prather et al.(1987), which as- 
sumes that within a model grid box a fraction of tracer 
mass is directly adverted from the convection base to the 
top of the convectively unstable column. ( ontinuity then 
requires an instantaneous subsidence to move tracer from 
the top down to the next lower box, throughout the col- 
umn. Thus total mass is conserved in the column while 
redistributing t racer. The diffusive scheme was based upon 
the convective adjustment used in the (K'M for convective 
parameterization of heat and moisture and in this sense, it 
is internally consistent with the dynamics. The diffusive 
scheme is different from that of the direct relocation, in 
that it tests adjacent, points in a column and mixes be- 
tween them if an instability exists. 1 his process is doin' 
from the' ground up so it is possible for tracer to slowly 


‘diffuse’ up through the column. The two schemes are 
fundamentally different and will have significantly differ- 
ent impacts upon the transport of tracers with different 
lifetimes. 

3. Results. 

Due to its relatively short life time, the radon is on av- 
erage mostly concentrated over the land area. Figure 1 
shows the mass mixing ratio of radon at the surface, aver- 
aged over the second month, on which can be seen many 
features associated with the large scale circulation. In the 
tropics, zonal transport produced by the trade winds is ev- 
ident as well as the impact of the tropical circulation due 
to the presence of the high pressure belt in the southern 
hemisphere. The persistence of such an anticyclonic circu- 
lation over the tropical ocean produces surface winds fol- 
lowing the west coast of Africa and south America causing 
a strong gradient on the edge of both continents. Average 
surface concentrations of radon obtained over most of the 
northern continents range from 100-300 pCi/m 3 which fall 
in the range of what is currently observed (Liu et al.,1984). 
However, spatial and temporal variability of radon is such 
that comparison between observation and model results 
over specific regions would necessitate further analysis (Fe- 
ichtcr, 1988; Jacob, 1990). 

At mid-and high-latitude in the winter season, the im- 
pact. of convective' mixing is weak no matter which scheme 
is used. Tropical regions are' however strongly affected 
by the choice of parameterization and show highly local- 
ized spatial and temporal variations. To illustrate the 
differences between the schemes, figures 2b and c show 
time snap-shots of the radon distributions for the two dif- 
ferent parameterizat ions. for 7 T West at model day 33, 
00:00 GMT. Figure 2a shows, for comparison, the result 
achieved when no convert ive mixing is applied. The ‘diffu- 
sive 1 and non-couvective mixing variants art; very similar 
with onlv a slight spread of tracer due to diffusion from the 
tracer plumes generated by synoptic scale ascent. 1 bis in- 
dicates that the ‘diffusive' technique spreads tracer slowly 
to high altitudes but has little power to inject tracer higher 
due to the relatively short lifetime of radon. The ‘reloca- 
tion* scheme however is strikingly different producing a 
high level tracer maxima in the tropics right up to the 
tropopause and lower stratosphere. In addition, without 
convective mixing the model tracer (due to the source at 
the ground) builds up strong gradients in the lowest lev- 
els (Figure 2a). producing numerical problems at land-sea 
boundaries and above continents. 

The mean residence time of lead in tin' stratosphere 
is considerably higher than in tin' troposphere, where the 
lead is washed out by precipitation. Previous studies have 
shown that the quantity of lead in the stratosphere is a 
significant, fraction, of the total (Moore. 19<3) and could 
correspond to 33Vf ol the global mass of lead in the at- 
mosphere (Lambert, 1982). Our results show that such a 
reservoir can be approached by using the relocation tech- 
nique 
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Figure 1 : Radon 222 surface mass mixing ratios, (pci/m 3 
at STP), for Radon variant ii) ’Diffusive’ convective mix- 
ing, Monthly averaged data for February. 

which gives a stratospheric mass fraction of 25% at quasi- 
steady state. The fraction given by the diffusive mixing 
scheme is of the order of 5-10% which is significantly lower. 

In spite of the slowness of the diffusive scheme, it still 
contributes to the transport of species into the higher tro- 
posphere. Budget analyses show in fact that close to 30% 
of the transport of radon to the tropopause could be at- 
tributed to this process. To allow the tracer to go further 
into the stratosphere, an additional type of parameteriza- 
tion could be needed in that case to represent the exchange 
of tracer that occurs across the tropopause. 

Part of the aim of these experiments was to determine 
the capabilities and limitations of the spectral technique. 
Initial model experiments reveal the model conserves mass 
well, giving a maximum global mass variation of 0.1 per- 
cent in 2 months. Punning at 132 resolution the radon and 
lead tracers together with the convective mixing schemes 
give rise to strong gradients locally, regionally and glob- 
ally. Despite the noisy nature of the convection and the 
permanent strong land-sea gradients which arise 1 the to- 
tal global magnitude of t he numerical errors remain small, 
less than 0.25 percent of the global mass being used to 
remove negative values in the tracer field. For Radon the 
combined impact of the ground source and the convective 
mixing is a relatively severe test on the use 1 of the spectral 
advert ive scheme. Consequently sharp gradients and low 
background values will be produced and numerical prob- 
lems be more likely. A surprising bonus of the convective 
mixing has been its ameliorating effect upon the need for 
hole-filling in tin 1 model. As the result, of active convec- 
tion, the gradients of |{n 222 are no longer so sharp and 
trapping in the lower model levels is reduced. 

A number of t racers in fut lire si udies will be like Hu 222 
in having a source restricted to over land. I his creates a 


tracer distribution with steep horizontal gradients at the 
continental boundaries. Such a distribution must be mod- 
eled by the advection scheme, maintaining the gradients, 
to be accurate. Figure 1 illustrates the spectral model’s 
description of Rn 222 near the surface, with strong gradi- 
ents in the horizontal, with the tracer mass predominantly 
‘clinging 1 to the land surfaces. Future tracers will also have 
detailed geographical source distributions, e.g. dependent 
on terrain characteristics. The impact of similar sources 
varying on very small model scales (sometimes below the 
models 3 — A Ax limit ) will be partially tested as a result 
of the ‘noisy’ nature of the convective mixing which to the 
upper troposphere is acting like a highly variable source 
function in both space and time. With these types of ad- 
ditional requirement for adverting chemical tracers, the 
spectral technique is one among many alternative schemes 
proposed. (JCMs using the spectral technique* at low reso- 
lution cannot resolve; st rong gradients. Other alternatives 
exist including semi Lagrangian and second-order moment 
techniques which can provide more 1 accurate advection but 
normally at the expense 1 of more computational effort and 
sometimes poorer mass conservat ion properties. 

The radon distributions exhibit, st rong transport of trac- 
ers into the upper troposphere particularly in the tropics 
both via synoptic and subgrid motions. The consequences 
of this transport are important in any study of tracers 
whose sourer is the* troposphere but. whose* sink lie's in 
the stratosphere. With this type of transport scenario 
in mind the radon ami lead tracers have been examined 
to view how the* model move's tracers to the 1 upper tro- 
posphere* and across the* tropopause into the lower strato- 
sphere. The impact of the* ‘relocation - convective mixing 
and large* synoptic ascent upon (lie* redistribution e>f radon 
from (he* boundary layer to the upper troposphere is con- 
siderable*. It could have also an important ef[e*ct, in the 
long t<Tm. upon the* cross equatorial transport ot long lived 
t racer. 
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Figure 2 : Radon 222 mass mixing ratio constant, lon- 
gitude cross-sections, longitude 1 71°\\'est. for variants: 
a) No convective 1 tracer mixing, b) Diffusive convec- 
tive tracer mixing, and c) direct-Relocat ion convective 
tracer mixing, lor Day 33, 00:00 GM I . (Scaled by 1 .0c ) 
Vertical axis shows data on model hybrid levels, where 
lj = f(p.<T). 

1. Conclusions 

A level 0 online tracer-model is now complete and test- 
ing indicates that mass conservation and negatives aie not 
significant problems but. that a better advectiou scheme 
may be needed if sharp gradients produced dynamically 
and chemically arc 1 to be maintained aitmatel\. However 
the spectral technique with sufficient resolution is still a 
powerful and viable alternat ive. Sensitivity studies reveal 


the importance of the parameterization of convective mix- 
ing on tropospheric tracer modeling. The model’s convec- 
tive mixing has maximum impact over the equatorial con- 
tinents, with very small effect in the winter hemisphere 
higher latitudes. In terms of global impact, the param- 
eterization of convective mixing based on the relocation 
technique appears to be more efficient in representing ad- 
equately the transport of long lived tracers into the lower 
stratosphere. However, the local impact of both schemes 
over specific regions still needs to be assessed with more 
analysis. It has been also noted that convective mixing 
acts to reduce spectral problems due to its diffusive effect 
on tracer gradients. Strong synoptic vertical motions pro- 
duce important transport of radon into the tropical up- 
per troposphere, greatly increasing the potential pool of 
tracer available for transport into the lower siratosphere. 
The convective mixing parameterization used are rela- 
tively simple and questions on the degree of entrainment 
throughout the column, fraction of mass relocated and 
the height of the ensemble average convection within a 
GCM grid box remain outstanding. Future modeling is 
intended using alternative and more complex schemes in 
order to assess and improve our ability to parameterize 
and interpret subgrid mixing within global models. 
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